The dynamics of coherence between coupled and decoupled one-dimensional 

quasicondensates 
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We reproduce the sub-exponential decoherence of one-dimensional quasicondensates observed in 
recent experiments. Counter-intuitively, the quasicondensates may decohere even when stongly 
coupled, if the temperature is large enough or the peak density is low enough to allow significant 
density fluctuations. We also propose an experiment to investigate the growth of coherence between 
two initially incoherent quasicondensates. We predict that the coherence will rise on a much slower 
timescale, and the final coherence again depends strongly on the density fluctuations. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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Questions of coherence are central to the study of Bose- 
Einstein condensates (BECs). These questions are inter- 
esting from the perspective of fundamental science, but 
also are crucial for the development of future applica- 
tions, possibly based on atom chips or optical lattices Ij. 
Typically, these devices create very elongated clouds, 
and can realise a quasicondensate regime at tempera- 
tures above a characteristic temperature that may be 
much less than the transition temperature Tc 0, H, 0, 0] ■ 
In this regime, global phase coherence is lost, and the 
phase fluctuates over a characteristic lengthscale which 
is less than the longitudinal length of the cloud. It re- 
mains an open question whether these quasicondensates 
could be exploited in practical applications, although ex- 
periments have attempted to perform interferometry in 
this regime [6|. An understanding of the coherence be- 
tween coupled and decoupled quasicondensates, as well 
as physical insight into their dynamics, is an important 
step towards their future use in interferometers. These 
dynamics are even more intruiging because quasiconden- 
sates are often also quantum one-dimensional systems, 
which sets important restrictions on their behaviour. 

In this letter, we study a recent experiment which 
realised one-dimensional quasicondensates on an atom 
chip. The one-dimensional regime occurs when ksT and 
UoUp are both less than hujr [1, 0] , where T is the tem- 
perature of the cloud, Up is its peak density, ujr is the 
radial (high) trapping frequency, and Uq = Anha/m, in 
which a is the s-wave scattering length and m is the mass 
of a single atom in the cloud. For the parameters in the 
experiment, UoUp = 2 x 10"^" < ftw^ = 3 x 10"^" J. The 
highest temperature considered in this paper is 200 nK, 
at which fc^T = 3 x lO"'^" J. However, when considering 
the dynamics of splitting, merging and weak links, it is 
important to question what we mean by one-dimensional. 
Clearly, during a splitting process into two separate one- 
dimensional clouds, the originally one-dimensional qua- 
sicondensate passes through a two-dimensional interme- 



diate stage. We address this ambiguity by perform- 
ing simulations in which the one-dimensional nature of 
the quasicondensate is not assumed. We find that this 
two-dimensional state provides a route for equilibration, 
which may be suppressed in one-dimensional clouds 
In fact, counter- intuitively, two quasicondensates may de- 
cohere despite being strongly coupled, if the clouds con- 
tain sufficiently large density fluctuations. This effect 
may prevent one from attaining the adiabatic limit when 
mergingquasicondensates in a double well or optical lat- 
tice This effect also implies that the details of 
splitting and merging protocols may have a large impact 
on the final state of the cloud. For example, a split- 
ting or merging protocol which involves a long-lived cou- 
pled state may cause significant decoherence between the 
two quasicondensates, which could be observed as heat- 
ing due to vortex production 13, U]. Once a dynamical 
equilibrium between the two coupled quasicondensates 
has been reached, we may then recreate the decoherence 
dynamics ij, Il3l . |14| in the experiment by abruptly sev- 
ering the link. In this way, we obtain the sub-exponential 
decay observed experimentally j3l and predicted theoret- 
ically [ll,[ll. We also propose an experiment to investi- 
gate the growth of coherence between two initially inco- 
herent quasicondensates when they are merged to create 
a link. At low temperatures, the coherence rises, but 
on a slow timescale, which we explain in terms of vor- 
tex dynamics, in an analagous effect to the Kibble-Zurek 
mechanism 15]. At higher temperatures, the presence of 
density fluctuations inhibits the growth of coherence. 

The results in this paper are obtained using the finite- 
temperature truncated Wigner method, which we devel- 
oped in a previous publication to analyse atom-chip in- 
terferometry of BECs [ll[ . This method has the advan- 
tage of capturing the zero-temperature mean-field de- 
pletion, which has a significant effect on high density 
elongated clouds. Essentially, the method models ther- 
mal fiuctuations by calculating the Bogoliubov excita- 



tions of the condensate mode, and then populating them 
with a thermal distribution (see Ref. [ll| for details). 
As in our previous paper we do not add excita- 

tions with an average occupation of less than one atom. 
The exclusion of vacuum fluctuations and hence inco- 
herent scattering allows us to clearly identify the phys- 
ical mechanisms of decoherence in our calculations, and 
also it eliminates any spurious decoherence caused by 
the numerical simulation of quantum noise in the trun- 
cated Wigner method To match with the exper- 
iments, we construct finite T initial states containing 
Nt = 5 X 10'^ ^^Rb atoms, and the trap frequencies in the 
axial and radial directions are cj^ = 27r x 5 rad s~^ and 
= 27r X 4000 rad s ^ respectively. Once we 
have prepared three-dimensional initial states, we evolve 
in time a two-dimensional slice in the y — Q plane, since 
the motion in the y-direction is frozen out. 

We begin by characterising a single cloud at K. Fig- 
ure HJa) shows the density profile in the y — Q plane 
of the central section (approximately one third) of the 
cloud. The figure illustrates the extreme aspect ratio 
of the system. Such clouds are commonly imaged af- 
ter expansion, during which they expand rapidly in the 
radial direction. Figure, [ijb) shows the full cloud af- 
ter 8 ms expansion. For the parameters in the experi- 
ment, = 15A^T (^^) V32Atfcs « 0.5 uK [3,@1, where 
jjL = 3.5 X 10~^° J is the chemical potential. This temper- 
ature is much less than Tc = 330 nK. Consequently, we 
would expect to find a true BEC with global phase co- 
herence at nK. To test this, we calculate the coherence 
of the cloud ^'i(r), where r is distance in the z-direction 
from some arbitrary point (0, 0, z) [axes inset in Fig.[l](a)] 
near the center of the cloud. This is done by comparing 
the phase 6{r + r„) to 0(r„), where r„ are N positions 
near the center of the cloud. We may define 5'i(r) math- 
ematically as 
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The resulting plot of \E'i(r) is shown in Fig. [He). The 
curve shows that coherence drops from one to almost zero 
over ^ 40 /im. This is much less than the longitudinal 
length of the cloud, which we calculate to be ~ 250 /im. 
We conclude that even the zero-T mean-field depletion 
is sufficient to destroy global phase coherence in this ex- 
tremely dense and elongated cloud. 

Although quasicondensates have no global phase co- 
herence, density fluctuations are suppressed due to the 
mean- field interaction 0, d, For example, the quasi- 
condensate shown in Fig. [UJa) has a reasonably smooth 
density profile. Density fluctuations may be observed fol- 
lowing expansion, when interactions become negligible, 
as shown in Fig. [Ub). However, as the temperature is 
increased, small density fluctuations may appear in-trap. 
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FIG. 1: (a) Typical atom density profile of central region 
(~ 1/3 of total cloud) of quasicondensate in y = plane 
(axes inset) at nK. (b) Corresponding atom density profile of 
complete cloud after 8 ms expansion. Upper [lower] horizontal 
bar indicates scale in (a) [(b)]. (c) Coherence 'i'l as a function 
of distance r in the z-direction within the quasicondensate at 
nK. (d) In-trap density fluctuations Fn as a function of T. 



We may characterise these density fluctuations with the 
quantity F„ (T), deflned as 



FniT) 



Jj.J^^n^z,T)dz- ^J^^niz,T)dz 



2L 



J^^n{z,T) dz 



(2) 

where n{z,T) is atom density at position (0,0,2;) and 
temperature T, and 2L is a distance less than the axial 
length of the cloud. We plot (T) in Fig. [BJc) for 2L = 
20 /im. The curve shows that density fluctuations are 
suppressed at low temperatures, but grow rapidly when 
the temperature is raised above 100 nK. 

We now simulate the dynamic equilibrium of two cou- 
pled quasicondensates [l8|, separated by a distance D, as 
a function of temperature. The coherence ^'2(t) between 
the two clouds as function of time t is measured as 



(3) 



where d (z) is the relative phase between the two clouds, 
and 2L is 100 /tm. This quantity is always averaged over 
five simulations with different initial conditions. We be- 



gin our simulations with two identical quasicondensates 
but with a very small amount of added noise in order t( 
provide a seed for the dynamic evolution. This noise is oi 
average much less than half a particle per mode, whicl 
was added in previous work to model incoherent scatter 
ing [l^, 0, 21 1 . Consequently, incoherent scattering i; 
prohibited in our calculations. We check this conditio! 
is satisfied by calculating ^'2(t) for D = 1.5 ^m, whicl 
is sufRciently large to suppress coupling between the twi 
clouds. In this case, the coherence does not drop signifi 
cantly from one over 80 ms, as shown in Fig.[2I^a) for 101 
nk (dashed curve) and 200 nk (dot-dashed curve). O 
course, in a real experiment incoherent scattering wouh 
cause decoherence, but our purpose here is to prepan 
coupled quasicondensates in dynamical equilibrium. Ti 
do this, we set _D to 1 /im, so that the height of thi 
barrier is roughly equal to the chemical potentials of thi 
clouds. Counter-intuitively, the coupling increases thi 
rate of decoherence. This effect is shown by the solic 
line in Fig. [2ja) at T = 100 nK. At this temperature 
the effect is relatively small, and the coherence plateau 
at a value of about 0.8. However, the effect is mucl 
more dramatic at 200 nK, where the coherence drops t( 
about 0.2 [dotted curve in Fig. [2l[a)]. It is intruiging tha 
this coupled system allows a route to equilibration, whicl 
may be suppressed for one-dimensional clouds [7] . Thesi 
results show that, suprisingly, there is no coherence fo 
coupled quasicondensates above a certain temperature 
and hence the experiments would have been impossibli 
at higher temperature, despite being well below Tc- 

The decoherence process occurs via small initial differ- 
ences in phase or density between the two clouds driving 
small Josephson oscillations, causing a transfer of atoms 
between the two wells. As atoms are gained or lost, the 
mean-field interaction in each quasicondensate responds 
by modifying the density profile. This leads to greater 
differences in phase or density between the two clouds, 
driving greater Josephson oscillations. However, this pro- 
cess does not ultimately lead to a complete loss of co- 
herence. Since the Josephson oscillations are incoherent 
over the phase coherence length of the quasicondensate, 
they create fluctuations in density over this distance. Ul- 
timately, the amplitude of the Josephson oscillations is 
restricted by the amplitude of the density fluctuations 
at the temperature of the cloud. At 100 nK, the den- 
sity fluctuations are small, as shown in the upper inset of 
Fig.[2l[a), which is an enlargement of the central region of 
the two coupled clouds at t = 50 ms. These small fluctua- 
tions should be compared to the much larger fluctuations 
in the corresponding image at 200 nK in the lower inset 
of Fig. [H^a). As previously shown in Fig. [TJc), the mag- 
nitude of the density fluctuations increases rapidly as the 
temperature is raised above 100 nK. 

Our calculations reveal the suprising result that sig- 
nificant decoherence may occur while the quasiconden- 
sates are strongly coupled, and consequently the details 
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FIG. 2: (a) Coherence ^2 plotted as a function of time for 
two coupled/uncoupled quasicondensates at 100 nK (solid 
curve/dashed curve) and 200 nK (dotted curve/dot-dashed 
curve). Upper (lower) inset shows the atom density profile in 
the central region of two coupled clouds in the y = G plane 
(axes inset) at t = 50 ms and T = 100 (200) nk. The field-of- 
view is 7.5 x 3 /xm. (b) Atom density profile in central region 
of a quasicondensate at 200 nK at t = 60 ms, after performing 
merging protocol A. The field-of-view is identical to that in 
(a) insets, (c) As (b), but for merging protocol B. 



of splitting and merging protocols may have a large im- 
pact on the cloud dynamics. To illustrate this, we con- 
sider two protocols in which two initially identical qua- 
sicondensates are merged from D = 1.0toZ? = 0.5 /im. 
In protocol A, the clouds are merged over 1 ms and then 
held for a further 59 ms. In this case, there is little exci- 
tation of the cloud, as shown in Fig.[2jb), which is an en- 
largement of the central region of the cloud at i = 60 ms. 
In protocol B, the quasicondensates are held dX D = 1.0 
/im for 59 ms, allowing them to reach dynamical equlib- 
rium, before being merged to D = 0.5 /im. In this case 
the density profile at t = 60 ms has large density fluc- 
tuations, and the enlargement of the central region of 
the cloud in Fig. ^c) contains two vortices (white cir- 
cles in center of figure). This figure shows that vortex 
production, as observed previously in atom interferome- 
ters [13, EH, can contribute to heating of the clouds, 
despite them being individually one-dimensional. 

Figure [2l[a) shows that our simulations of coupled qua- 
sicondensates reach a dynamical equilibrium after 80 ms 
evolution. For i ^ 80 ms, our simulations capture the 
correct variations in the occupations of the Bogoliubov 
modes of two coupled quasicondensates. We may then 
simulate the decoherence of the two quasicondensates by 
abruptly splitting them to D = 1.5 /im and consequently 
severing the link between them, allowing the Bogoliubov 
modes in each cloud to evolve freely. Hence we avoid the 
challenging problem of directly modelling the uncertain- 



FIG. 3: (a) Decoherence of two quasicondensates at 100 nK 
after splitting to -D = 1.5 (dashed curve) and 1.15 (solid 
curve) ^m. (b) Dashed curve: double logarithmic plot of 
dashed curve in (a) showing sub-exponential decay of coher- 
ence. Solid line: line of best fit. 



ties in the splitting process in the experiment The 
resulting loss of coherence at 100 nK is shown by the 
dashed curve in Fig. [31 which has a very similar timescale 
and form to that reported in experiment. We also carry 
out a partial splitting to £) = 1.15 /um, and observe the 
coherence saturate at ^ 0.6 (solid curve). 

Theory has predicted that the coherence should decay 
as e(-*/*o)°, where a = 2/3 [H, [13]. We investigate 
this by plotting In (— In (^'2)) against In (t), as shown in 
Fig. [Sfb). The gradient of the line of best fit (sohd hne) 
is 0.7, in rough agreement with theory and experiment. 
We also extract a value of 23 ms for to from the same 
graph, in good agreement with the value predicted by 
the analytic results in Ref. [l^. 

We now extend our analysis to coupling of two ini- 
tially incoherent quasicondensates. We find that when 
the two quasicondensates are held at D = 1.0 fim at 100 
nk, the coherence rises over about 150 ms (solid curve 
in Fig. [4]). Crucially, the growth of coherence when two 
quasicondensates are coupled occurs on a wildly differ- 
ent timescale to the loss of coherence when two quasi- 
condensates are decoupled (shown in Fig. [3]). This can 
be understood in terms of the underlying physical pro- 
cesses. When two quasicondensates are decoupled, the 
unequal splitting of the thermally occupied Bogoliubov 
modes causes decoherence on a relatively fast timescale. 
When two quasicondensates are coupled, vortices form 
in low density region between the two clouds at posi- 

3 El. 



tions where the relative phase happens to be tt 
These vortices are pinned by the barrier potential, so co- 
herence may only occur once they have either migrated 
in the z-direction to the edges of the cloud, or collided 
with vortices of opposite rotation and annihilated. This 
is relatively slow process, occuring over a timescale of 
150 ms for the parameters of the experiment. This 
effect is analgous to the Kibble-Zurek mechanism [l5| . 
in which vortices may survive for long timescales after 
a rapid quench across a phase transition. In contrast to 
the results at 100 nK, at 200 nK the coherence remains 



FIG. 4: Establishment of coherence when two initially inco- 
herent quasicondensates are coupled at 100 nk (solid curve) 
and 200 nk (dotted curve). 



low (dotted curve in Fig. |4|). This is again due to the 
density fluctuations which occur above 100 nK. 

In summary, we have examined the coherence of cou- 
pled and decoupled quasicondensates, characterised the 
time-dependent dynamics and identified the underlying 
physical processes. It is important to note that, although 
we focus on the coherence of two quasicondensates in a 
double well, our findings can also be applied to splitting 
and merging in optical lattices. Crucially, our results 
show that there will be heating when quasicondensates 
are merged in a two-dimensional optical lattice jprevent- 
ing one from attaining the adiabatic limit [3, Q , if the 
temperature is high enough or the density low enough to 
allow significant density fluctuations. 
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